Linear response theory around a localized impurity in the pseudogap regime of an 
anisotropic superconductor: precursor pairing vs the d-density-wave scenario 
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We derive the polarizability of an electron system in (i) the superconducting phase, with 
d-wave symmetry, (ii) the pseudogap regime, within the precursor pairing scenario, and (iii) the 
d-density-wave (dDW) state, characterized by a d-wave hidden order parameter, but no pairing. 
Such a calculation is motivated by the recent proposals that imaging the effects of an isolated 
impurity may distinguish between precursor pairing and dDW order in the pseudogap regime of 
the high- T c superconductors. In all three cases, the wave- vector dependence of the polarizability is 
characterized by an azymuthal modulation, consistent with the d-wave symmetry of the underlying 
state. However, only the dDW result shows the fingerprints of nesting, with nesting wave-vector 
Q — (n, 7r), albeit imperfect, due to a nonzero value of the hopping ratio t' jt in the band dispersion 
relation. As a consequence of nesting, the presence of hole pockets is also exhibited by the (q, w) 
dependence of the retarded polarizability. 
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I. INTRODUCTION 

Imaging of the electronic properties around an isolated 
nonmagnetic impurity such as Zn in the high-T c super- 
conductors (HTS) has provided direct evidence of the 
unconventional nature of the superconducting state in 
the cuprates, and in particular of the d-wave symmetry 
of its order parameter below the critical temperature T c 
[111 H B In the underdoped regime of the HTS, vari- 
ous models have been proposed to describe the pseudogap 
state above T c . 

Several experimental results provide substantial evi- 
dence of a pseudogap opening at the Fermi level in un- 
derdoped cuprates for T c < T < T* , even though no 
unique definition of the characteristic temperature T* is 
possible, as it generally depends on the actual experi- 
mental technique employed (see Refs. H, E3 for a review, 
and refs. therein). Also, the doping dependence of T* 
is still a matter of controversy [jj. Owing to its d-wave 
symmetry, the pseudogap has been naturally interpreted 
in terms of precursor superconducting pairing. In partic- 
ular, the pseudogap has been associated to phase fluctu- 
ations of the order parameter above T c @ (see Ref. @ 
for a review). Within this precursor pairing scenario, 
the phase diagram of the HTS can be described as a 
crossover from Bose-Einstein condensation (in the un- 
derdoped regime) to BCS sup erconductivity (in the over- 
doped regime) 
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Recently, it has been proposed that many properties of 
the pseudogap regime may be explained within the frame- 
work of the so-called d-wave-density scenario (dDW) 
[l3L IT3 , Il5j . This is based on the idea that the pseudogap 
regime be characterized by a fully developed order pa- 
rameter, at variance with the precursor pairing scenario, 
where a fluctuating order parameter is postulated. The 
dDW state is an ordered state of unconventional kind, 
and is usually associated with staggered orbital currents 
in the Cu0 2 square lattice of the HTS [H [H EH . 
Much attention has been recently devoted to show the 
consistency of the dDW scenario with several experimen- 
tal properties of the HTS EH- These include transport 
properties, such as the electrical and thermal conductiv- 
ities I20l l2l| an d the Hall effect HHH, thermodynamic 
properties |24l l25j . time symmetry breaking |2q|. and 
angular resolved photoemission spectroscopy (ARPES) 
|27| . The possible occurrence of a dDW state in micro- 
scopic models of correlated electrons has been checked in 
ladder networks psj . 

It has been recently proposed that direct imaging of the 
local density of states (LDOS) around an isolated impu- 
rity by means of scanning tunneling microscopy (STM) 
could help understanding the nature of the 'normal' state 
in the pseudogap regime The idea that 

an anisotropic superconducting gap should give rise to 
directly observable spatial features in the tunneling con- 
ductance near an impurity was suggested by Byers et al. 
[331, whereas earlier studies [34| had considered pertur- 
bations of the order parameter to occur within a distance 
of the order of the coherence length £ around an impu- 
rity. Later, it was shown that an isolated impurity in 
a d-wave superconductor produces virtual bound states 
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close to the Fermi level, in the nearly unitary limit |35| . 
Such a quasi-bound state should appear as a pronounced 
peak near the Fermi level in the LDOS at the impurity 
site [3^. as is indeed observed in Bi-2212 and YBCO 

S3 

In the normal state, the frequency-dependent LDOS 
at the nearest and next-nearest neighbor sites, with re- 
spect to the impurity site, should contain fingerprints of 
whether the pseudogap regime is characterized by pre- 
cursor pairing (3^1 or dDW order [23, E| • This is due to 
the fact that while pairing above T c without phase co- 
herence is a precursor of Cooper pairing, and therefore 
of spontaneous breaking of U(l) gauge invariance, the 
dDW state can be thought as being characterized by the 
spontaneous breaking of particlc-holc symmetry, in the 
same way as a charge density wave breaks pseudospin 
SU(2) symmetry [33]. The LDOS around a nonmagnetic 
impurity in both the dSC, the dDW and the competing 
dSC+dDW phases in the undcrdoped regime has been 
actually calculated e.g. by Zhu et al. j^j. 

In this context, a complementary information is that 
provided by the polarizability F R (q,uu) of the system, 
which gives a measure of the linear response of the charge 
density to an impurity potential. In the case of d- 
wave superconductors, it has been demonstrated that the 
anisotropic dependence of the superconducting order pa- 
rameter on the wave-vector q gives rise to a clover-like 
azymuthal modulation of F R (q,u>) along the Fermi line 
for a 2D system [io|. 

These patterns in the q dependence of F R (q,oj) are 
here confirmed also for a more realistic band for the 
cuprates. In addition to that, the dDW result also shows 
fingerprints of the Q = (-7T, ir) nesting properties of such 
a state. 

The paper is organized as follows. In Sec. [01 we review 
the expression of the polarizability for the d- wave super- 
conducting state (dSC) and derive that of the <i-wave 
pseudogap regime, within the precursor pairing scenario 
(dPG). In Sec. IIIII we derive the polarizability for the 
dDW state. By allowing nonzero values of the hopping 
ratio t' jt in the dispersion relation [ill l42l E^ j. we will 
explicitly consider the case in which perfect nesting is de- 
stroyed. Such a case is relevant for the study of the dDW 
state, given its particle-hole character. In Sec. IIVI we 
consider the competition of dDW order with an subdom- 
inant dSC state in the underdoping regime. In Sec. [V] 
we present our numerical results for the polarizability in 
the dSC, dPG, and dDW states, both in the static limit 
and as a function of frequency. We eventually summarize 
and make some concluding remarks in Sec. IVII 



II. LINEAR RESPONSE FUNCTION IN THE 
dSC AND dPG STATES 

Within linear response theory, the displaced charge 
density 6p(r) by a scattering potential V(r) in the Born 



approximation is given by 

Sp(r) = J V(v')F K (r-T',E F )dv', (1) 

which implicitly defines the linear response function 
F R (r, Ep) at the Fermi energy E-p. Here and in the fol- 
lowing we set the elementary charge e = 1. Its relevance 
in establishing the electronic structure of isolated impuri- 
ties in normal metals and alloys has been earlier empha- 
sized by Stoddart et al. 0, EH ■ In momentum space, 
Eq. Q readily translates into 6p(q) = V (q) F R (q, E F ) , 
showing that, for a highly localized scattering potential 
in real space [V(r) = Vo<5(r), say], the Fourier transform 
<5p(q) of the displaced charge is simply proportional to 
F R (q,£ F ). 

In the presence of superconducting pairing, the gener- 
alization of the linear response function is given by the 
density-density correlation function (polarizability) [4^ | : 

F(q,z^) = Tri^l^r 3 a(k, iW „) 

" u„ k 

x T 3 G(k - q, iw n - iujy) (2) 

where {?(k, ioj n ) is the matrix Green's function in Nambu 
notation, (3 = T _1 is the inverse temperature, uj v = IvitT 
is a bosonic Matsubara frequency, r, are the Pauli matri- 
ces in spinor space, the summations are performed over 
the N wave-vectors k of the first Brillouin zone (1BZ) and 
all fermionic Matsubara frequencies w„ = (2n + l)7rT, 
and the trace is over the spin indices. Here and be- 
low we shall use units such that K = ks = 1 and lat- 
tice spacing a = 1. The retarded polarizability is de- 
fined as usual in terms of the analytic continuation as 
F R (q, oj) = F(^ioj l , 1 ► lj ~\~ i0 + ). In the normal state, 
Equation J5J correctly reduces to the Lindhard function 
for the polarizability of a free electron gas |47j . 

In the following, by specifying the functional form of 
Q in the case of pairing with and without phase coher- 
ence, we will in turn derive the explicit expression for 
F in the superconducting phase with a o?-wave order pa- 
rameter, and in the pseudogap regime, characterized by 
fluctuating d-wave order (precursor pairing scenario). 

A. Superconducting phase 

We assume the following BCS-like Hamiltonian: 
H d sc = X]£ kC ks c i« + v ^' c U c -ki c -^'l c ^'h ( 3 ) 

ks kk' 

where c^ s (ck s ) is a creation (annihilation) operator for 
an electron in the state with wave-vector k and spin 
projection s G {T> 1} along a specified direction, and 
£k = et — /i, with ek the single-particle dispersion re- 
lation: 

ek = — 2i(cos k x + cos k y ) + At' cos k x cos k y , (4) 
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where t = 0.3 cV, t'/t — 0.3 are tight binding hopping 
parameters appropriate for the cuprate superconductors, 
and n the chemical potential. In Eq. © , Hck' is a model 
potential, which we assume to be separable and attractive 
in the d x 2_ y 2-w&ve channel: Vkk' = — A^kSk', with = 
i(cosfc x — cos k y ) and A > 0. Under these assumptions, 
the Hamiltonian, Eq. @, is characterized by a nonzero 
superconducting order parameter (cktc-ki), leading to 
a nonzero d-wave mean-field gap Ak = A gk below the 
critical temperature T c . 

Making use of the explicit expression for the matrix 
Green's function GdSC m the superconducting state [isj : 



r n ■ \ iuj nTo + CkT 3 + A k n 
5dSc(k ' ia; " ) = {^^El ' (5) 

with Ek = (££ + A 2 .) 1 / 2 the upper branch of the super- 
conducting spectrum and tq the identity matrix in spin 
space, and performing the trace over spin indices and the 
summation over the internal frequency |47j in Eq. (j2J), we 
obtain the linear response function for a d-wave super- 
conducting state |46j : 



-Pdsc(q, iv v ) 



N ^ 



(likUk-q - Wk«k-q 



2 //(^k)-/(i?k- q ) 



-(tik-qVk + Mk^k-q) 



Ek — Ek-q — iu v 
2 //Cgk)+/(^k-q)-l 



H.c, 



V ^k 



H.c. 



k-q 



(6) 



where u 2 . = |[1 + £k/-Ek], = 1 — u 2 . are the usual 
coherence factors of BCS theory, and /(e) = (1 + e' 36 ) -1 
is the Fermi function at temperature T. In the limit of 
zero external frequency and T — > 0, Eq. © reduces to 
the static polarizability studied in Ref. |4jj for a d-wave 
superconductor. 

B. Pseudogap regime, within the precursor pairing 
scenario 

In the pseudogap regime, for T c < T < T* . within the 
precursor pairing scenario |9(, one assumes the existence 
of Cooper pairs characterized by a 'binding energy' Ak 
having the same symmetry of the true superconducting 
gap below T c , but no phase coherence. In other words, 
no true off-diagonal long range order develops, and one 
rather speaks of a 'fluctuating' order ||. This means 
that the quasiparticle spectrum E^ = + A 2 .) 1 / 2 is 
still characterized by a pseudogap Ak = A <7k opening at 

I 



where we are implicitly assuming T c < T < T* . 



the Fermi energy with ci- wave symmetry, but now with- 
out phase coherence. Therefore, the diagonal elements of 
the matrix Green's function GdPG coincide with those of 
its superconducting counterpart, Eq. (JHJ, while the off- 
diagonal, anomalous elements are null: 

The effects due to a finite lifetime of the precursor Cooper 
pairs can be mimicked by adding a finite imaginary 
energy lincwidth iV to the dispersion relation entering 
Eq. (JJJ, or by substituting the spectral functions as- 
sociated with the quasiparticle states with 'broadened' 
ones, as discussed in Appendix 1X1 The relation between 
the two approaches and with analytical continuation has 
been discussed in Appendix A of Ref. . 

Within this precursor pairing scenario, Equation J2J) in 
the pseudogap regime then reads: 



(8) 

I 

III. LINEAR RESPONSE FUNCTION IN THE 
dDW STATE 

The mean-field Hamiltonian for the d-density-wave 
state is |l3j : 

#dDW = X]K kC ks Cks + iD kcl s Ck + Q s ], (9) 
ks 



N 



,2„,2 



( f(E k ) - /(£ k -q) 
01 \Ek- £k-q - iOJ v 



H.c. 



2 2,22 
k-q«k + UuV 



k^k-qj 



/(£ k ) + /(gk- 

El- + E\r-CI - 



- 1 



H.c. 
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where the summation is here restricted to all wave- 
vectors k belonging to the first Brillouin zone, Q = (tt, 7r) 
is the dDW ordering wave-vector, and Dy = D g^ is the 
dDW order parameter. As anticipated above, the dDW 
state is characterized by a broken symmetry and a well- 
developed order parameter, at variance with the precur- 
sor pairing scenario of the pscudogap regime. Such a 
state is associated to staggered orbital currents circulat- 
ing with alternating sense in the neighboring plaqucttes 
of the underlying square lattice. As a result, the unit 
cell in real space is doubled, and the Brillouin zone is 
correspondingly halved. At variance with other 'density 
waves', the dDW order is characterized not by charge or 
spin modulations, but rather by current modulations. 

The nonzero, singlet order parameter $q breaks pseu- 
dospin invariancc in the particle-hole space: 



J 



Cks') = i$QgkS s 



(10) 



Whereas it possesses c?-wave symmetry, as expected, its 
imaginary value leads to the breaking of a relatively 
large number of symmetries, such as time reversal, par- 
ity, translation by a lattice spacing, and rotation by 7r/2, 
although the product of any two of these is preserved (see 
Ref. [21j for a detailed analysis). 



Introducing the spinor tyl = (c^ 



Hamiltonian. 
MM: 



), the dDW 

Eq. ©, can be conveniently rewritten as 



ks U k+Qs 



#dDW = 2 *LK e k - -") r + £k T 3 + AcTi^ks, (11) 
ks 

where = ^(ek ± ek+cj), and the prime restricts the 
summation over wave- vectors k belonging to the reduced 
('magnetic') Brillouin zone only. Notice that £k+Q = 
=FeJ. Correspondingly, the matrix Green's function at 
the imaginary time r can be defined as <?dD\v(k, r) = 
-(T r * ks (r)^j cs (0)}, whose inverse reads III, IP : 



^dDW 



(k, iu n ) 



lUJ n 



6, 



6c+Q 



(12) 



In the case of perfect nesting (t' = 0) for the dispersion 
relation, Eq. Q), Sharapov et al. [2lJ explicitly find 

g dDW {k,iuj n ) = — — — — — — 3 — ^ — i ( 13 ) 



(iuj n + nY 



to be compared and contrasted with Eq. (JSJ) for the su- 
perconducting phase. Notice, in particular, the different 
way the chemical potential [i enters the two expressions. 

In the general case (t' ^ 0), perfect nesting is lost, and 
we have to refer to the general form of C^w, Eq. (|14f) . 
One finds: 



£dDw(k, iu) n ) 



1 



(iw n - E+)(iuj n - £ k ) 

— £k+Q 



k 



(14) 



where E^r = — fi+e^± y (e k ) 2 + D 2 . are the two branches 
of the quasiparticle spectrum obtained by diagonalizing 
Eq. © HI- Notice that E^ +Q = E±. 

In the limit t 1 = 0, Eq. I|14|) correctly reduces to 
Eq. H13|) . even though it is not straightforward to ex- 
press Eq. (|14f> in the same compact matrix notation. 
In the limit of perfect nesting (t' = 0), the disper- 
sion relation, Eq. i s antisymmetric with respect to 
particle-hole conjugation, £k+Q = — £k- As a result, 
E^ = —\i ± (e k + D 2 .) 1 ^ 2 , which is to be contrasted with 
the quasiparticle spectrum of the superconducting state 
or the pseudogap state within the precursor pairing sce- 
nario, ±_Bk = ±[(ek — M) 2 + ^k] 1 ^ 2 - The difference comes 
again from the fact that the Bogoliubov excitations in 
the dSC and the dPG states are Cooper pairs, while the 
dDW ordered state is characterized by particle-hole mix- 
ture 0112. 

The polarizability in the dDW state is derived in Ap- 
pendix^ Wc just quote here the final result, which can 
be cast in compact matrix notation as: 

^dDw(q,^) = Tr 4 i^yi «£ dDW (k, iuj n ) 

" Un k 

x K<?dDw(k - q, iujn ~ iw v ), (15) 

where now k = tq + t± , and G<idw is the matrix Green's 
function for the dDW state, Eq. i|14f> . Performing the 
frequency summation |47l |. one eventually finds: 



1 



E 



/(££)-/(££-,) 



N k i^EL-K-*-™* 



(16) 



IV. COMPETITION BETWEEN dSC AND dDW 
ORDERS 

In the underdoped regime, it has been predicted on 
phenomenological grounds that the dDW order should 
compete with a subdominant dSC phase [l3j . This has 
been confirmed by model calculations at the mean-field 
level |25l 1^3, showing that indeed an existing broken 
symmetry of dDW kind at high temperature suppresses 
that critical temperature for the subdominant dSC or- 
dered phase. Recently, the competition between dDW 
and dSC orders has been shown to be in agreement with 
the unusual T-dcpcndcncc of the restricted optical sum 
rule, as observed in the underdoped HTS |5l| . 

In order to take into account for the competition be- 
tween the dSC and dDW orders at finite temperature, 
one has to separately consider the electron states within 
the two inequivalent halves of the Brillouin zone. There- 
fore, it is convenient to make use of the 4-components 
Nambu spinor ^ = ^Tkj,)i or explicitly: 



* k = 



( CkT \ 

C k+QT 

VcLk_ Qi / 



(17) 
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At the mean-field level, the Hamiltonian for the compet- 
ing dSC and dDW phases thus reads 



dSC+dDW 



(18) 



where 7J k is the 4x4 Hermitean matrix defined by 



H k = 



( a iD k A k 
iDk £ k +Q -A k 
M -£_ k »D_k 



(19) 



V -A* -iD_ k -£_ k _q 
with real eigenvalues Eki (here, i = 0, . . . 3) given by 
E li = A^ + £>k - £k£k+Q 

+ (£k + £k+Q 



Ek) 2 



(20) 



and orthonormal eigenvectors |ki) (for given k in the re- 
duced 1BZ). It may be straightforwardly checked that 
Eki reduces to the dSC superconducting spectrum, ±Ek, 
and to the dDW quasiparticle dispersion relations, E^, 
in the limits D D = (pure dSC) and A = (pure dDW), 
respectively, when the halving of the 1BZ is removed. 

In the particle- hole symmetric case (t 1 = 0), mak- 
ing use of the nesting properties described in Sec. IIIII 
Eq. I2U|) simplifies as: 



£k* 



±1 



Dl 



(21) 



For [i — 0, the four branches of the spectrum degenerate 
into the Dirac cone I52ll53l 



-Ek. 0=2, 1=3 — i 



A k 



(22) 



thus showing that the two gaps have the same role, i.e. 
the system may be equivalently described a s a dDW or 
a dSC superconductor, with a d-wavc gap y/ A 2 . + in 
either case. Either a non-zero hopping ratio (t' /t ^ 0) or 
a hole-doping away from half-filling (fj, ^ when t' = 0) 
destroys this particular symmetry, and one has to resort 
to the eigenvalues Eki, Eq. (j20J). 

In order to obtain the Green's functions in the 
dSC+dDW case, it is useful to introduce the 4x4 ma- 
trices (i, j = 0, ... 3) 



whose algebra is given by 



(23) 



(24) 



Then the matrix Hamiltonian Eq. (|19f) takes the form 

H k = e k r^3 + e k+Q r33 - ^r 30 - L> k r 02 + A k ri 3 , (26) 

whence the inverse Green's function (now a 4 x 4 matrix 
in Nambu space) straightforwardly follows as 

£dSC+dDw( k ' iu} n) = «w„r 00 - Hk- (27) 

As in the dDW case, Eq. l[T5|) (see also Appendix IB")) , the 
linear response function in the dSC+dDW case can be 
given a compact matrix form as 

fdsc+dDw(q,«<v) = Tr - tp ££dsc+dpw(k, iu n ) 

" u n k 

x Kf/ dsc +dDw(k - q, iuJ n - iu v ), (28) 

where now the vertex matrix in the 4x4 Nambu spinor 
space is given by 



k = T3 <x> k = 



r a + T l 

-To - T\ 



(29) 



Finally, it can be shown that Eq. (|28|l also admits the 
following spectral decomposition, analogous to Eq. (|TT)|) : 



-FdSC+dDw(q, i&v) 
1 

N 



k i,j 



Eki — Ek-qj — lLO u 



(30) 



where Uy(k, q) = Tr (£-P k iK-P k -qj), and Pki = |ki)(ki| is 
the orthonormal projector operator on the i eigenstate of 
the matrix Hamiltonian, Eq. (|19|) . 



V. NUMERICAL RESULTS AND DISCUSSION 

We have evaluated numerically the polarizability for 
the dPG and the dDW cases, Eqs. © and JTJJ, and for 
the mixed dSC+dDW case, Eq. I|30|) . as a function of the 
relevant variables. Our numerical results for the pure 
dSC case turn out to be very similar to the dPG case (at 
least over the range of variables considered below), and 
will not be shown here. In the dPG and in the pure dDW 
cases, we adopt the following set of parameters, which are 
believed to be relevant for the cuprate superconductors: 
t = 0.3 cV, t' /t = 0.3, [i = —t, corresponding to a hole- 
like Fermi line and a hole doping ~ 14.3%, A Q = D Q = 
0.06 eV in the dPG and in the dDW cases, respectively 
E3, andT= 100 K. 



where e^u is the totally antisymmetric Levi-Civita ten- 
sor, and 



(25) 



A. Zero external frequency 

In order to make contact with earlier work |40j , we first 
consider the case of zero external (bosonic) frequency, 
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(a) (b) 



FIG. 1: Static polarizability for the dPG state in momen- 
tum space, F d pG(q, 0), Eq. JSJ [in eV -1 , panel (a)], for 
A = 0.06 eV, T = 100 K, and fj, = -t = -0.3 eV, cor- 
responding to the hole-like Fermi line in the normal state 
shown in panel (b). Panel (b) also reports the special points 
F = (0,0), X = (it, 0), M = (it, it), along with the points A 
and B where the Fermi line £k = intersects the symmetry 
contour T-X-M-Y. 




-2 2 K 



(a) (b) 

FIG. 2: Static polarizability for the pure dDW state in mo- 
mentum space, -FdDw(c[, 0), Eq. 11611 [in eV -1 , panel (a)], 
for Do = 0.06 eV, T = 100 K, and fi = -t = -0.3 eV, 
corresponding to the hole pockets E^ — shown in panel 

(b) . Panel (b) also reports the special points T — (0,0), 
X = (7r,0), M — (vr,7r), along with the points Bi and Bi 
where the line E^ — intersects the symmetry contour V— 
X-M-V. 



lu u = 0, in the timc-ordcrcd polarizabilitics, Eqs. JHJ and 

Our numerical results for the wave-vector dependence 
of F(q, 0) over the 1BZ in the dPG and in the pure dDW 
cases are shown in Figs. ^andj^ respectively. As a result 
of the d-wave symmetry of both the pseudogap within the 
precursor pairing scenario, and of the dDW order param- 
eter, F(q, 0) is characterized by a four-lobed pattern or 
azymuthal modulation |4jj. However, the dDW case is 
also characterized by the presence of 'hole pockets', cen- 



tered around Q/2 = (ir/2,ir/2) and symmetry related 
points, due to the (albeit imperfect) nesting properties 
of the dDW state, with nesting vector Q = (it, it) (sec 
Fig. [2b) ■ Such a feature is reflected in the q dependence 
of i 7 dD\v(qi 0), which is characterized by local maxima at 
the hole pockets, for the value \i = — t of the chemical 
potential considered here. (Other values of the chemi- 
cal potentials give rise to analogous features, which are 
absent in the dPG case.) 

Fig. [3] shows our numerical results for the static polar- 
izability F(q,0) in the mixed dSC+dDW case, Eq. (50} . 
Representative values of the amplitudes of the dSC and 
dDW order parameters have been taken as in Ref. [Hcj . 
viz. A = O.lt = 0.03 eV, D Q = 0.08t = 0.024 eV, 
at T = O.Olt ~ 35 K, for a particle-like Fermi line 
in the underdoped regime (i = 0.3 eV, t'/t = 0.3, 
/x = -0.2016 eV). Panel (b) of Fig. shows the con- 
tour plot of the eigenvalue spectrum E^, Eq. (|2*U|l . The 
latter is characterized by pronounced minima near the 
hot spots at Q/2, which evolve into cone-like nodes in 
the limit of pure dSC (D = 0), or in the very special 
case t' = (i = [see Eq. (|2*2*)) ]. Accordingly, Fig. 0a, for 
the static polarizability -F(q, 0) over the 1BZ is charac- 
terized by local maxima at the hot spots centered around 
Q/2, as is the case in the pure dDW case (cf. Fig- Et)- 
Whereas the precise behavior of the static polarizability 
F(q, 0) is of course determined by the actual amount of 
dSC+dDW mixing at a given temperature and doping, 
we can conclude that a sizeable dDW component man- 
ifests itself through the appearance of hole pockets cen- 
tered around Q/2 in the q-dependence of F(q, 0), also in 
the presence of a dSC condensate. 

We have next evaluated the spatial dependence of 
F(r,0) (not shown), by Fourier transforming F(q, 0) to 
real space. While F(r, 0) is characterized by Friedel- 
like oscillations as |r| increases from the impurity site, 
as expected 0i l50j . these radial, damped oscillations 
are superimposed by an azymuthal modulation, due to 
the d-wave symmetry of the normal state, both in the 
dPG and in the dDW cases. As a consequence, F(r, 0) is 
characterized by a checkerboard pattern, closely related 
to the symmetry of the underlying square lattice, with 
local maxima on the nearest neighbor and local minima 
on the next- nearest neighbor sites. Since these features 
are common to both the dPG and dDW cases, the spa- 
tial dependence of the charge density oscillations is not 
directly helpful in distinguishing between the dPG and 
dDW states. However, real-space and wave-vector de- 
pendences of several quantities of interest for STM stud- 
ies can be easily connected by means of Fourier trans- 
form scanning tunneling microscopy (FT-STM) tech- 
niques (see, e.g., Ref. |54j . and refs. therein). Such 
a technique has been proved very effective in detecting 
large-amplitude Friedel oscillations of the electron den- 
sity on the Be(0001) [HE^ and Be(1010) surfaces H3, 
and has been recently discussed in connection with ex- 
perimental probes of fluctuating stripes in the HTS [58j . 
In particular, FT-STM experiments |55l l56l l57j have evi- 
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(a) (b) 

FIG. 3: Static polarizability for the mixed dSC+dDW state in 
momentum space, -FdSC+dDw(q, 0), Eq. I30H [in eV~\ panel 
(a)], for A = O.lt = 0.03 eV, D = 0.08t = 0.024 eV, T = 
O.Olt ~ 35 K l50ll . and /j, — —0.2016 eV, corresponding to a 
particle-like Fermi line closed around the F point (underdoped 
regime). Panel (b) shows the contour plots of the eigenvalue 
spectrum E ki , Eg. 1201 

denced the role of correlation and reduced dimensionality 
in establishing such 'giant' Friedel oscillations in the elec- 
tron density. 

B. Frequency dependence 

We have next evaluated the frequency dependence of 
the retarded polarizabilities, in both the dPG and the 
dDW cases. Figures 0] and [S] show our numerical results 
for the uj dependence of Fjp G (q, uj) and F r dw (c[,uj), re- 
spectively. Each curve refers to either the real or the 
imaginary part of -F R (q, w) as a function of uj, for a fixed 
value of wave-vector q along the symmetry contour T- 
X-M-T in the 1BZ (see Figs.[IJ) and E for its definition). 
While q runs along such contour, the Fermi line £k = 
is traversed twice (once at point A along X-M and once 
at B along M-T, for the hole-like Fermi line considered 
here; see Fig. ^0, while the hole-pocket contour defined 
by — is traversed twice along the M-T line (points 
Bi = B and B 2 in Fig. Eh). 

As a consequence of the summation over k in either 
the full or reduced BZ in Eqs. © and (|16|l . respectively, 
ReF R (q,w) is an even function of uj, while ImF R (q, w) 
is an odd function of u> in both the dPG and dDW cases. 
Therefore, the different contributions of particle and hole 
states in the two cases is averaged out, and no asymmet- 
ric peaks in the u> dependence of such quantities are to 
be expected in the dDW case, as is the case for the local 
density of states [^j, ■ 

On the other hand, the existence of hole pockets cen- 
tered around Q/2 in the dDW state is clearly responsible 
for the different u dependence of RcF^q (Fig.0^,) versus 
R° ^dDW (Fig- 130 j sa y ; as Q runs along the T-X-M-T 
contours. While Re F^ PG is characterized by a single rela- 



tive maximum for uj > for all wave- vectors q under con- 
sideration, Re -F^w possesses two relative maxima (or a 
relative maximum and a shoulder) for u> > 0. These two 
maxima tend to merge into a single one for Bi -< q -< B% , 
i.e. inside the hole pocket, and for q m A, i.e. at the 
intersection of the free-particle Fermi line with the X-M 
side (Fig. . Likewise, the single relative maximum for 
uj > in Rc-F^^r shifts towards larger frequencies as q 
runs from T to X, is 'diffracted' at A along the Fermi line 
as q runs from X to M, and 'bounces back' at B, again 
along the Fermi line, as q runs from M back to T. A 
similar analysis may be performed for ImF R in the two 
cases (Figs.^J) and|5p). 

As in the static limit, the competition of a sizeable 
dDW order parameter with an underlying dSC conden- 
sate does not give rise to qualitatively different results 
in the w-dependence of the polarizability, with respect to 
the pure dDW case. 

One may conclude that, in both the dPG and dDW 
cases, the evolution with q of the features in uj depen- 
dence of F R (q,uj) are closely related to the location of 
wave-vector q with respect to the Fermi line, and may 
therefore serve to indicate the presence of hole-pockcts, 
as is the case for the dDW state. 



VI. CONCLUSIONS 

Motivated by recent STM experiments around a local- 
ized impurity in the HTS, we have derived the polariz- 
ability (density-density correlation function) -F(q, uj) for 
the pseudogap phase, both in the precursor pairing sce- 
nario and in the d-density-wave scenario. Expressions for 
the same function have been derived also in the under- 
doped regime, characterized by competing dSC+dDW 
orders. 

In the static limit (here defined as the limit of zero 
external frequency for the time-ordered correlation func- 
tion), the q dependence of F(q, 0) reflects the d-wave 
symmetry of the precursor pairing 'pseudogap' or of the 
dDW order parameter, with an azymuthal modulation 
consistent with a clover-like pattern, as expected also for 
a superconductor with an isotropic band |40j . However, 
at variance to the dPG case, the q dependence of the 
static polarizability in the dDW state clearly exhibits the 
presence of hole pockets, due to the (albeit imperfect) 
nesting properties of the dDW state, with nesting vec- 
tor Q = (tt,7t). Qualitatively similar results to the pure 
dDW case are obtained also in the mixed dSC+dDW, 
thus showing that hole pockets are a distinctive feature 
of dDW order. Such a behavior is confirmed by the r 
dependence of the static polarizability in real space. A 
detailed comparison with experimental data for the in- 
dependence of the charge density displacement would of 
course require a much more detailed knowledge of the q 
dependence of the impurity potential, here crudely ap- 
proximated with an s-wavc Dirac (5-function. In particu- 
lar, the presence of higher momentum harmonics in the 
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FIG. 4: Frequency dependence of the real [panel (a)] and 
imaginary parts [panel (b)] of the retarded polarizability, 
-F R (q, lj), in the dPG case, for wave- vector q varying along a 
symmetry contour T-X-M-T in the 1BZ (see Fig. QJ>) . All 
curves have been shifted vertically for clarity, by an amount 
proportional to the path length from F to actual wave-vector 
q along such symmetry contour (see right scale). Dotted line 
is the zero axis for F R (0, u>). All other parameters as in Fig.Q 



impurity potential may break the d-wave symmetry of the 
possible correlated or ordered states (dPG, dSC, dDW) 
here studied. Also, an extension of the present Born ap- 
proximation for the impurity perturbation, e.g. to the 
T-matrix formalism, would afford a more reliable com- 
parison with experimental results. 

An analysis of the frequency dependence of the re- 
tarded polarizability F R (q,u>) reveals that the q evolu- 
tion of the features (local maxima or shoulders) in the 
lu dependence of this function is closely connected with 
the relative position of wave- vector q with respect to the 
Fermi line, and is therefore sensitive to the possible pres- 
ence of hole pockets, as is the case for the dDW state. 
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FIG. 5: Same as Fig. 3] but for the dDW case. Special points 
in the 1BZ are as in Fig. 
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APPENDIX A: FINITE LIFETIME EFFECTS 

In order to take into account for finite lifetime effects 
on the linear response function for the pseudogap regime 
within the precursor pairing scenario, we write the diag- 
onal elements of the matrix Green's function as 

fe(k,z^») = ^ / du; Au{Ku} \ (Al) 

Z7T J lU) n — LU 

where Au(k,u) = 2Tru^S(Ek — uj), A 2 2(k, lo) = 
2irv%.5(Ek +lo) are the appropriate spectral functions for 
BCS theory. 

A finite energy lincwidth T can be attached to the en- 
ergy state by replacing the ^-functions in the spectral 
functions An with broader ones, e.g. a Lorcntzian func- 
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tion a(u) = ^T/(uj 2 + T 2 ). Setting 

A(k,uj) =2%[uia(E k -uj)+vla(E k + u)], (A2) 
in the static limit one obtains: 

^dPG*.^ w °xt - uj - - / — - — — — 

J IjJ — LO 

x </>(q, lj, uj') dw du', (A3) 

where 

<f>(q,w,w') = -^i^4(k )W )yl(k-q, W '). (A4) 

APPENDIX B: POLARIZ ABILITY FOR THE 
dDW STATE 

In order to derive the analog of the polarizability, 
Eq. |PJ, for the dDW state, we start with considering 
the density-density correlation function: 

F(q,T) = -(T r p(q,T)p(-q,0)>, (Bl) 

where p(q, r) = J2k s c q S ( r ) c k+qs(T) is the electron den- 
sity operator, and T T denotes ordering with respect to 
the imaginary time r. Application of Wick's theorem 
then yields 

F(q,r) = ^(T rCk+qs (r)4, s ,(0))(T TCk ,_ qs ,(0)cL(r)) 

kk' 

-(/7(q,0))(p(-q,0)), (B2) 



the last term being a constant with respect to r, which 
can be neglected in Fourier transforming to the Matsub- 
ara frequency domain. In the dDW state, the contribu- 
tions of terms like Eq. I|1U|I should be explicitly consid- 
ered. Therefore, we make use of the identity 

E-^EVk + Zk+q), (B3) 

k k 

for the summations on both k and k' in Eq. i|B2|) , where 
the prime restricts the summation to wave-vectors k be- 
longing to the reduced (magnetic) Brillouin zone. Af- 
ter Fourier transforming to the Matsubara frequency do- 
main, one eventually has: 

-FdDw(q,«^) = 4y^^y^ [GniKiVn) + £12 (k, iu n ) 
' w n k 

+g 2 \ (k, iuj n ) + 5 2 2(k,io;„)] 

x [<7n(k - q, iw n - iuj v ) + £12 (k - q, iu n - iu u ) 

+Q21 (k - q, iu> n - iujy) + 5 2 2(k - q, iuj n - iu v )\, 

(B4) 

where Qij are the entries of 5dD\v(k, iuJ n ) m Eq. (|14|) . 
The last expression can then be cast into the compact 
matrix form, Eq. I|15f) . by introducing the constant aux- 
iliary matrix k = tq + t\ = |^ . Equation l|B4(l sim- 
plifies further, by observing that C/12 = — £/2i and that 
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